Formation of van der Waals molecules in buffer gas cooled magnetic traps 



N. Brahms,!' 2 rj. y_ Tscherbul,^' 3 p. Zhang,^ J. Klos,^ H. R. 

Sadeghpour,^ A. Dalgarno,^'^ J. M. Doyle, ^'^ and T. G. Walker*^ 

^Department of Physics, University of California, Berkeley, California 97420 
^Harvard-MIT Center for Ultracold Atoms, Cambridge, Massachusetts 02138 
^ITAMP, Harvard- Smiths Oman Center for Astrophysics, Cambridge, Massachusetts 02138 
Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742 
^Department of Physics, Harvard University, Cambridge, Massachusetts 02138 
^Department of Physics, University of Wisconsin-Madison, Madison, Wisconsin 53715 

(Dated: April 29, 2011) 

We show that a large class of helium-containing cold polar molecules form readily in a cryogenic 
buffer gas, achieving densities as high as 10^^ cm""^. We explore the spin relaxation of these 
molecules in buffer gas loaded magnetic traps, and identify a loss mechanism based on Landau- 
Zener transitions arising from the anisotropic hyperfine interaction. Our results show that the 
recently observed strong thermal dependence of spin change in buffer gas trapped silver (Ag) is 
accounted for by the formation and spin change of Ag^He, thus providing evidence for molecular 
formation in a buffer gas trap. 



Techniques to produce, trap, and manipulate cold neu- 
tral and ionic molecules hold great promise for new dis- 
coveries within chemistry and physics. Cold molecules of- 
fer new ways to improve precision spectroscopy and mea- 
surement IJ, and explore quantum collective phenomena 
such as quantum magnetism [2\ and high-Tc supercon- 
ductivity 3J. The long-range anisotropic interactions of 
polar molecules in optical lattices might be tuned with 
electric fields, allowing for the design of robust quantum 
information processing [4j and quantum simulations of 
condensed-matter systems [5]. 

To date, two types of cold trapped molecules have 
been produced in the laboratory. The first are molecules 
formed via photo- and magneto-association from ultra- 
cold atomic gases, creating ultracold alkali dimers [7] 
and trimers |8]. The second type are stable dipolar 
molecules slowed or cooled from high temperatures using 
Stark deceleration W or buffer gas cooling [TD] . Here we 
introduce a third class of trappable molecules — van der 
Waals (vdW) complexes — and provide compelling theo- 
retical support for the formation of Ag^^He heteronuclear 
molecules in a recent experiment, in which Ag atoms were 
magnetically trapped in the presence of a dense '^He gas 

Weakly bound vdW complexes play a key role in chem- 
ical reaction dynamics, surface interactions, and non- 
linear optical phenomena in dense atomic and molecu- 
lar vapors [I2HI6] . In particular, vdW complexes con- 
taining He atoms have attracted special attention due to 
the existence of quantum halos in ''He2 [II] and simi- 
lar exotic many-body states in ^Hes [H]. Due to their 
small binding energies and vulnerability to collisions. He- 
containing complexes have previously only been observed 
in the collision-free, transient conditions of supersonic 
jets [H]. 

Here we show how a wide variety of cold, trappable 
vdW molecules can be produced in buffer gas loaded 



magnetic traps. We provide a general model for the 
formation of vdW pairs and calculate expected molec- 
ular densities for a sampling of magnetically trappable 
species. Using ah initio calculations of interaction 
potentials and hyperfine interactions, quantum cal- 
culations of collision dynamics, and detailed Monte 
Carlo simulations, we show that recent spin relaxation 
measurements [11 provide convincing evidence for the 
formation of AgHe molecules. 

Molecular formation kinetics. We begin by outlining the 
general model for the formation of A"He molecules via 
three-body collisions of a species X (atoms, molecules, 
or ions) with He atoms: 

X + He + He ^ AHe + He, (1) 

where K and D are the rate constants for three-body 
recombination and collision-induced dissociation. From 
this equation, we obtain the formation and dissociation 
kinetics: 

"XHc = -nx = nx/n - rixHc/Td, (2) 

where Ui denotes the density of species i, and the forma- 
tion time Tf and dissociation time obey 1/rf = Kn^^ 
and l/rd = D riHe ■ 

In thermal equilibrium, nxHc = implies nxHc = 
K(^)nx?^Hc, where the chemical equilibrium coefficient 
k(T) = K/D can be evaluated from statistical thermo- 
dynamics [T7]: 

- = — - iv^) E 9.e-''^ ^ (3) 

where /cb is Boltzmann's constant, T is the temperature, 
/J, is the reduced mass of the AHe molecule, and — is the 
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TABLE I: Predicted ground-state energies eo and molecular 
population ratios "^"^ = KriHe of a few species compatible 
with buffer gas loaded magnetic traps |18H20) . 

X^He X*He 
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State 


— eo 


"XHo b 
"x 


— eo 


"XHo i> 
n-x 
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''•S'3/2 


2.13 
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•5*3/2 


2.70 


84 


3.42 


18 


Cu 


^•5*1/2 
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0.015 


1.26 


0.008 


Ag 


•5*1/2 


1.40 


0.16 


1.85 


0.06 


An 


^•5l/2 


4.91 


3 X 10^ 


5.87 


1 X 10^ 


"Energies 


in cm^^, 


1 cm~^ R 


s 1.4 K. is the dissociation energy. 



''Population ratios are evaluated for njje = 3 X 10^® cm ^, at 
300 inK for X^He and at 700 mK for X^Re molecules, using only 
the level with energy eg- 

dissociation energy of molecular state i, with degeneracy 
gi. Eqn. ([s]) shows that rixHo increases exponentially as 
the temperature is lowered, or as the binding energy of 
the molecule is increased. A large number of vdW com- 
plexes have binding energies that are comparable to or 
greater than 1 K, and thus formation is thermodynami- 
cally favored in buffer gas cooling experiments. Table [l] 
gives a sampling of candidate species for molecular pro- 
duction and predicted population ratios "^^^ — ktiho at 
standard temperatures and buffer gas densities. 

The molecular density will come into equilibrium on a 
timescale Tc, determined by rf and t^- From Eqn. ([2]): 

= + Tf-i = D riHc (1 + KnHo). (4) 

An estimate for the rate of approach to chemical 
equilibrium is obtained by assuming that D ^ trwe'^/'^^"^, 
with a a typical low-temperature gas kinetic cross 
section and the exponential factor giving the relative 
probability that a colliding He atom has enough en- 
ergy to dissociate the molecule. Choosing worst-case 
values cr = 10" cm^, v = 50 m/s, T = 0.3 K, and 
riHe = 10^^ cm^'^ gives Tc < 20 ms for species with 
I Co I < 1.4 cm~^, and less than 400 ms for all eq. 
Typical buffer gas trap lifetimes of seconds, therefore, 
indicate that the molecule population reaches equilib- 
rium even for those species having |e|/fc_B larger than 1 K. 

Spin stability. We investigate the spin stability of 
magnetically trapped molecules, considering those 
formed from S'-state atoms for simplicity. The Hamilto- 
nian for the ground electronic and vibrational state of 
an XHe(I]) molecule in a magnetic field of strength B is 

^moi = Cat + Axlx-S+B -(2^38 + flxlx + MHoIho) + 

2 

7N • S + AhcIhc • S + cW— ^ n%(^)[lHc ® S](2), 

q=-2 

(5) 



where e^r are the rotational energy levels, N is the ro- 
tational angular momentum, S is the electron spin, 7 is 
the spin-rotation constant, Ix is the nuclear spin of X 
with moment ^x, -^x is the atomic hyperfine constant, 
and fj,B is the Bohr magneton. The last two terms in 
Eqn. ([5]) describe the isotropic and anisotropic hyperfine 
interaction of S with a ■^He nuclear spin Ihc- We neglect 
the weak anisotropic component of the interaction. 

The typical trap used for buffer gas cooling is a mag- 
netic quadrupole, with \B\ — at the center, linearly 
increasing to a few Teslas at the edge. The weak-field- 
seeking states of both the atoms and the molecules are 
stably confined in such a trap. 

Collisions with free He atoms can cause loss of the 
trapped XHe molecules 21j. The coherence of the pre- 
cessing electronic spin is lost in the collision, and the 
internal N-S and Iho^S interactions can cause spin re- 
laxation, with the Zeeman energy released into the ki- 
netic energy of the collision complex. However, the B 
field strongly decouples the electron spin from other in- 
ternal angular momenta. The probability to suffer spin 
relaxation in a collision due to an interaction of the form 
aJ • S, after averaging over the trap distribution of mag- 
netic fields, is p = ^^24jk]^T^- Estimates of 7, Aho, and 
c, all on the order of a few MHz [22], give p ^ 10~^, far 
too small to be a significant factor for trap lifetimes. 

We have identified a spin-relaxation mechanism, due 
to adiabatic transitions at avoided level crossings in- 
duced by the anisotropic hyperfine interaction, explained 
in depth in our discussion of AgHe below. While this 
is the dominant loss process for certain X^He pairs, we 
find that trap lifetimes > 1 s are still achievable for 
S-state vdW heteronuclear molecules formed in buffer 
gas traps. 

AffHe molecules. We will now sketch our analysis 
for the case of Ag'^He. Recently, we reported on the 
trapping of Ag atoms in a "^He buffer gas [iT]. The 
observed Ag spin-change rate decreased by almost two 
orders of magnitude as temperature was increased from 
300 to 700 mK; the data are reproduced in Fig. [l] 

To analyze the trap loss of Ag'^He, we first calculated 
the potential energy curve to an expected accuracy of 
10%, shown in Fig. [2] (a). We employed the partially- 
spin restricted coupled cluster method with single, dou- 
ble and perturbative triple excitations (RCCSD(T)) ^23] 
as implemented in the MOLPRO suite of programs [24]. 
The potential supports one vibrational bound state, with 
iV = 0,1,2 at energies = -1.40,-1.04,-0.37 cm-\ 
and a quasibound N = 2, state at £3 — 0.48 cm^^, with a 
calculated lifetime of 1 ns. The molecular hyperfine con- 
stants were evaluated as a function of internuclear dis- 
tance r using a highly correlated density functional the- 
ory approach |25j . The values, averaged over the ground 
state molecular wavefunction, are = —0.9 MHz and 
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FIG. 1: Experimentally observed Ag-'^He spin change rate co- 
efficients (squares) , compared with the atomic spin-relaxation 
model (dashed line), and molecular spin change rate coef- 
ficients, obtained from Monte Carlo simulations (solid line, 
shaded area indicates 68% C.I.). The simulations were used 
to fit for the Ag^He binding energy at the 380 mK datum, 
yielding eo = — 1.53±0.08 cm^^. Also shown is a comparison 
to the chemical equilibrium coefficient k{T) (dotted line, right 
axis). 



c — 1.04 MHz. The Zeeman spectrum of AgHe, consist- 
ing of identical hyperfine manifolds for each rotational 
level, is shown in Fig. [5] (b). 

We estimate the Ag'^He formation rate constant K by 
assuming that pairs first form in the N — S scattering 
resonance, and are subsequently quenched by rotational 
relaxation in additional collisions with He. Using the 
resonant recombination model described in [26], we cal- 
culate values for K between 0.8 and 1.0 x 10^"^^ cm^/s, 
for T between 0.3 and 0.7 K, giving an equilibration time 
of < 10 ms for all T and rifjo in the experiment. 

To calculate the spin relaxation rate due to Ag-'^He 
collisions, we performed scattering calculations over the 
experimental temperatures and magnetic fields using the 
rigorous quantum formalism developed in |251 127) . We 
set an upper bound on the spin-rotation constant 7 < 
1.6 MHz by scaling up the general formula from ^22} to 
match the observed relaxation at 0.7 K. We then averaged 
the relaxation rate coefficients over the trap field and 
thermal collision energy distribution [281 ES] to produce 
the dashed curve shown in Fig. [T] 

The calculated rate constant shows only a weak tem- 
perature dependence and exhibits marked disagreement 
with the experimental data below 0.6 K, suggesting 
that the observed loss rates cannot be explained by 
atom-atom collisions. We have also considered the 
Ag^He-^He coUisional relaxation mechanism and, as 
expected from the estimates above, find it far too small 
to explain the experimental results (see |25|). 



Adiabatic transitions We finally consider the possi- 
bility of adiabatic transitions resulting from avoided 
crossings between the low-field-seeking = and 
high-field-seeking N = 2 Zeeman levels, shown in 
Fig. [2](c). At a magnetic field of 1.06 T, the anisotropic 
hyperfine interaction couples these levels, causing 
avoided crossings ,30J with splittings ranging from 155 
to 380 kHz. Trapped AgHe molecules crossing this 
magnetic field can experience a spin flip, causing the 
molecule to be expelled from the magnetic trap. We 
model this process as a Landau-Zener type adiabatic 
transition with probability 

Mip = 1 - exp — , (6) 

which depends on the coupling strength V , the molecu- 
lar translational velocity and the magnetic field gradi- 
ent V-B. plz versus B — \v ■ VB\ is shown in Fig. [2| (c). 

To understand how these transitions lead to an anoma- 
lous loss of Ag atoms, we use a direct simulation Monte 
Carlo approach [551 ISI]- For each temperature, 40,000 
trapped Ag and Ag^Heare simulated under the effects 
of the magnetic trapping field and hard-sphere collisions 
with the ^He gas, with elastic collision probabilities cal- 
culated using the ab initio AgHe potential energy in 
Fig. [2] (a). Each collision has a probability to cause 
molecular formation or dissociation according to the rate 
coefficients K and D = K/k, respectively. We fit the re- 
sults of the simulation to the measured spin change rate 
at 380 mK, extracting a value for cq of — 1.53±0.08 cm~^, 
indicating a theoretical underestimation of the poten- 
tial depth on the order of 10%. Using this fit potential, 
we then simulated the spin change rate vs. temperature, 
finding agreement to the data over two orders of magni- 
tude without further fitting. 

We note that under typical experimental conditions, 
n-x, nxHc ^ ^iHc, and the formation rate of X2 dimers is 
negligible. Moreover, the experimental temperature far 
exceeds the 1 mK binding energy of '*He2, ruling out the 
possibility for the formation of He dimers, and tiho is too 
low to produce the large He clusters that are found in 
supersonic expansions [35]. For many species ArHep>i 
clusters may also form. The equilibrium density of these 
clusters is set by thermodynamics and can be made neg- 
ligible by raising the system temperature. 

In summary, our kinetic model for vdW pair forma- 
tion and relaxation shows that weakly bound atom-He 
heteronuclear molecules can form in copious quantities 
in buffer gas cooling experiments at sub-Kelvin temper- 
atures. By performing quantum calculations of Ag-'^He 
interactions and dynamics, and identifying the trap loss 
mechanism for vdW molecules, we have explained the 
marked, heretofore unexplained, thermal dependence of 
buffer gas trapped Ag spin relaxation. 
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FIG. 2; (a) Ab initio Potential energy for Ag^He (^E+), showing the three bound A*' = 0, 1, 2 molecular states (solid lines) and 
the quasibound N — 3 state (dashed line), (b) Zeeman level spectrum of Ag'^He (circle indicates region of avoided crossings), 
(c) Avoided level crossing between the | A, mjv, ms, mjAg, m/He) = |0, 0, 1/2, —1/2, 1/2) and |2, 2, —1/2, —1/2, —1/2) states. 
The inset shows the spin flip probability as a function of B for V = 380 kHz crossings (solid line) and 190 kHz crossings 
(dashed line). 



The molecules we predict will form in buffer gas ex- 
periments are trappable, and are slightly polar, with 
dipole moments ~ 0.02 D. It may be possible to pre- 
vent dissociation of the molecules by quickly reducing 
the He density — cryogenic valves can remove He in 40 ms 
[33]. Were this performed with 300 mK Ag at a density 
of 10""^^ cm~'^, chemical equilibrium would be lifted at 
7iHc < 10^^ cm^^, yielding a trapped sample of Ag'^He 
at a density of 7 x 10^" cm"'^. 
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